Scaling in a simple model for surface growth in a random medium 
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a: 

Surface growth in random media is usually governed by both the surface tension and the random local 
forces. Simulations on lattices mimic the former by imposing a maximum gradient m on the surface heights, 
and the latter by site-dependent random growth probabilities. Here we consider the limit m — > oo, where the 
, surface grows at the site with minimal random number, independent of its neighbors. The resulting height 
distribution obeys a simple scaling law, which is destroyed when local surface tension is included. Our model is 
equivalent to Yee's simplification of the Bak-Sneppen model for the extinction of biological species, where the 
<—j . height represents the number of times a biological species is exchanged. 
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Invasion percolation Q is a model for viscous fingering, in which at each time step a segment 
. of the interface moves into the capillary channel which imposes the (overall) minimum resistance. 



Tj" | Numerically, this is modeled by a lattice of sites, where each site is assigned a random number between 
■ zero and one, and by moving the interface into the perimeter site with the smallest random number. 
. After a while, most of the perimeter random numbers are relatively large (above the percolation 
CN | threshold) , and growth occurs in "bursts" which explore the vicinity of where they started j| . 

Another family of growth models is based on the "solid on solid" concept, where one ignores 
"£3 ■ overhangs and allows only steps which increase the height h of the interface (relative to the initial line 
£h or plane). Such models, which also allow for surface tension, include the Kardar-Parisi-Zhang (KPZ) 
' j [3 1 model, where the random resistance depends on horizontal position x and on time, and the Parisi 
[I] model, where this resistance is local, depending on both x and h. A discrete version of the latter, 
O ■ by Kim and Kosterlitz [g], only allows steps which obey the constraint \h(x) + 1 — h{x ± 1)| < m, 
J-^ | with m = 1. Sneppen (and Jensen) Q considered variants of this model, in which neighbors continue 
to be updated until the Kim-Kosterlitz constraint is obeyed everywhere. Bak and Sneppen j7|] then 
^ . used a related model to describe the extinction of biological species. In their model, they always 
^ | updated also the nearest neighbors of the growing site, irrespective of the height gradients. However, 
unlike Sneppen and Jensen, they did not continue these updates in an "avalanche" beyond the nearest 
neighbors. The "bursts" were thus more localized than in the original Sneppen model. 

Recently, Yee ^ introduced a simplified version of the Bak-Sneppen model, which can be identified 
as the m — > oo limit. In this model, there is no surface tension constraint on the growth, and thus each 
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"pillar" continues to grow as long as the random number in front of it is smaller than those facing all 
other pillars. Although this model preserves the long range information on all the perimeter "resis- 
tances" , it omits the interactions between neighbors, and is thus geometry and dimension-independent. 
In the language of invasion percolation, this model could describe a collection of independent capillary 
channels, penetrated from one end by a viscous fluid from the same reservoir, and with independent 
randomly varying resistances along each channel. The removal of interactions allowed an analytic 
solution of many aspects of this model ||, in agreement with simulations, which were phrased in the 
biological language. In the present paper we re-interpret this model as one of surface growth, and 
test if the resulting height distribution of this growth obeys a standard scaling law. We start with the 
definition of the model and its simulation and end with a partial analytical treatment. 

The Yee version of the Bak-Sneppen model takes an array of L random numbers rj, i = 1, 2, . . . , L 
initially distributed between zero and unity. At each time step t — > t + 1, the smallest of the L random 
numbers is replaced by a new random number. After some time ||, nearly all random numbers are very 
close to unity, and growth occurs in large "bursts" where a single "pillar" grows, similar to invasion 
percolation [§]. This version of the model f|| is independent of any geometry; the full Bak-Sneppen 
model replaces also the lattice neighbors of the lowest random number, and thus depends on the 
assumed lattice geometry; it gives a threshold x c below unity such that after a long time x c < < 1 
for nearly all i. The random numbers rj can be interpreted as fitness j?]] in biological evolution, or as 
the quality of court decisions in a judicial system with law-by-precedent ||. We now use the physical 
interpretation Q. 

With every element i we associate a height variable hi which initially is zero for all i, and then 
increases by unity every time this element i gets a new random value (which corresponds to the 
new perimeter site in front of the moving interface). We can imagine a deposition process in which 
bricks drop down onto the site with the lowest random number rj and then change this rj into another 
random number, or as the viscous fluid moving along the i'th (one dimensional) capillary channel, 
and reaching a new resistance. What is the probability distribution function (proportional to the 
histogram) of the observed heights? For this purpose we stop the growth whenever the highest hi 
value reaches a predetermined value L z (as usual in invasion percolation simulations and experiments 

I)- 

To avoid overcrowding of our figures we binned the observed heights into powers of two; that means 
the A:-th bin contained heights between 2 k ~ 1 and 2 k — 1. It is plausible that for 1 <C hi < L z and large 
L z the results should depend mainly on the ratio hi/L z . Figure 1 shows in its three parts the binned 
histogram N(h) in the scaled form N/L z versus h/L z , giving a good data collapse for large enough 
L; the three parts correspond to L z /L = 0.1, 1 and 10 and for large L seem to give the same curve. 
The initial linear increase, for h <C L, in these log-log plots, together with the fact that the bin size 
increases as the height h, would imply that without binning this increase corresponds to a constant 
probability distribution function oc N. At large heights we see a cut-off since h > L z is impossible. In 
fact, Fig. 2 shows that the unbinned distribution follows roughly an exponential, cx exp(— 6h/L). 

The time r after which the tallest pillar hits the top, hi = L, is shown in Fig. 3 for our squares 
(and rectangles). It increases roughly as I 2 , z ~ 1.8. However, a slight curvature suggests that 
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asymptotically the exponent z may be 2. Indeed, a plot of t/L 2 versus 1/L ' 3 (Fig. 4) seems to 
approach a finite limit for L — > oo. This means that when the tallest pillar hits the top, a finite 
fraction of the whole L x L lattice is occupied by the bricks of the pillars (or by the invading fluid). 
About the same intercept 0.07 was also found from the high and the flat rectangles of Fig. 1, for 
t/(L z L) versus 1/L a3 . 

As a function of time t, the average height H trivially always increases linearly in time, while the 
width W oc t@ of the height distribution has an exponent (3 increasing towards unity for t increasing 
towards r (not shown). 

Now an interaction between neighboring sites i is introduced. If i is updated, then also i + 1 
and i — 1 are updated if \hi — hi+i\ or \hi — respectively, are > m. Here m is a fixed number 

between and L z . This interaction corresponds to some sort of surface tension, which tries to avoid 
too large gradients in the height profile hi. It also makes our model one-dimensional, since now the 
neighborhood introduces a geometry. The limit m = corresponds to the Bak-Sneppen model (always 
updating of neighbors) and the limit m = L z to the simplified Yee version (no updating of neighbors). 

Figure 5 shows how the time r, the average height H =< hi >j and the surface roughness < 
(hi — H) 2 > x / 2 depend on this new parameter m; they go neither to infinity nor to zero, but the height 
H has a pronounced minimum at small m. Because of this new length m, the above simple scaling in 
terms of h/L z no longer works, even if as in Fig. 6 we take m to be that value (5 to 20) for which H 
has a minimum. Thus the non-interacting version obeys simple scaling while the interacting version 
depends on the geometry (here: one-dimensional only) and disobeys simple scaling. 

As stated, many features of the Yee model were calculated analytically by Newman |J. Specifically, 
at time t the smallest random number on the "perimeter" was shown to grow as 

approaching unity at long times. Indeed, if a "pillar" starts growing at a time t ^> L then it will 
continue to grow for n consecutive steps, during which the new random numbers encountered by this 
"pillar" are smaller than x(t). The probability of such an n-step growth was found to be exponential: 

Pn = X n ~ 1 (l-x) , (2) 

yielding an average step of length 

/ \ 1 t+L ^ 

{n) = —x = — ■ (3) 

Indeed, our simulations confirm "bursts" whose length grows linearly with time. Also, the same 
formalism yields 

<(n-(n)) 2 ) = ^±^ or ((„ - (n)) 2 )/(n) 2 = ^ (4) 

which also agrees with our simulations for large L and t, see Fig. 7. This unifractal distribution, in 
which n scales (for 1 -C L -C t) as n ~ t/L, is clearly different from that expected in other growth 
models mentioned in our introduction. 



3 



Note also that a given growth stops when the next random number is larger than x, which happens 
with probability 1—x = L/(L+t). Assuming that the random numbers are distributed equally between 
zero and one, every site will be encountered at least once when 1 — x becomes smaller than 1/L, i. e. 
at times of order 1? . It is not clear yet if this result relates to our numerical values for r, which were 
asymptotically consistent with being oc L x L z oc I? . 

Ref. H also derived the probability to find a growth "run" of length n at any time, 

OO i 

p n = wtr-Hi-sc*)) 2 *- . (5) 
U n 

The final height hi of a "pillar" is a sum over such "run" lengths, which grow longer and longer with 
time. Our simulations show that for L z = L, the typical number of such runs is of order 3-5. Thus, 
the distribution N(h) discussed above should in principle be a convolution of distributions like N n . We 
evaluated several such convolutions, with a variable number of "runs", and they all seem to converge 
asymptotically (for large h) to N(h) ~ 1/h, which differs from the exponential form found numerically 
in Fig. 2. At the moment we have no explanation for this discrepancy. 

We thank P.M.C. de Oliveira for drawing our attention on the Sneppen model, the German-Israeli 
Foundation for supporting our collaboration and the supercomputer center in Jiilich for time on their 
Cray-T3E. 
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Scaled histogram of pillar heights, Lz = L; L = 30 (diam.), 100 (+), 300 (sq.), 1000 (x), 3000 (tr.), 10000 (* 
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Figure 1: Scaling of the pillar height histograms: Curves for different L collapse, except for very 
small and very large heights. The shape of the L x L z lattice is a square (part a), a flat rectangle 
(part b), and a high rectangle (part c). Up to 640,000 samples were averaged over. 



5 



Sum over 50000 samples for unbinned unsealed height distribution, L = Lz = 1 00 and 300 
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Figure 2: Examples of unbinned and unsealed distributions of heights, showing a roughly exponential 
decay, for L = L z = 100 and 300. 
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Number of time steps to reach height Lz for length L = Lz (diam.), L=1 0*Lz (+), L=Lz/1 (sq.) 
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Figure 3: Time to reach the top, versus L, for squares, flat, and high rectangles. 
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Correction to simple scaling for time to hit the top, L = Lz = 30 to 1 0000 
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Figure 4: The times of Fig. 3 are shown to follow t/L = 0.07 + 0.5/L + . . . for large systems. 
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Time for highest pillar to hit the top; 1 000 samples 1 00 x 1 00; interacting pillars 
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Figure 5: Influence of interaction parameter m at L = L z = 100, averaged over 1000 samples, with 
m = L corresponding to the simplified Yee model and m = to the one-dimensional Bak-Sneppen 
model. Part a shows the time to reach the top, part b the height (diamonds) and the width (+) of 
the surface defined by the pillar tops. 
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Scaled histogram of interacting pillars, L = 30 (diam.), 100 (+), 300 (sq.), 1000 (+), 3000 (tr.); 1000 runs 
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Figure 6: Failed scaling of pillar height histograms for thousand L x L samples. The interaction 
parameter m was taken to give a minimum of H(m) (diamonds in Fig. 5b) and varies from 5 to 20. 
In contrast to the analogous Fig. 1, the different curves do not collapse to one curve. 
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(<n*n>-<nxn>)/(<n><n>) for L = Lz = 30000 (+, 10 samples) and 100000 (diam., 1 sample) 
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Figure 7: Comparison of < (n— < n >) 2 > / < n > 2 with theoretical time dependence t/(t + L), 
where n is the length of a stretch ("run" Q) of uninterrupted updatings of the same site. 



11 



